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Abstract. Hole and spin Green's functions of the two-dimensional t-J 
model of the Cu02 planes are calculated in an approximation which re- 
tains the rotation symmetry in the paramagnetic state and has no presumed 
magnetic ordering. In this approximation Green's functions are represented 
by continued fractions which are interrupted by decouplings with a vertex 
correction determined from the constraint of zero site magnetization in the 
paramagnetic state. Obtained results are shown to be in good agreement 
with Monte Carlo and exact diagonalization data. The calculated spectra 
demonstrate a pseudogap and a band splitting which are similar to those 
observed in Bi-based cuprate perovskites. 



For the description of strongly correlated electrons in Cu02 planes of 
perovskite high-T c superconductors the two-dimensional t-J model is widely 
used (for a review see Ref. [1]). One of the most fruitful analytical meth- 
ods for this model uses the spin-wave description of magnetic excitations. 
Some spectral peculiarities obtained in this approach are similar to those 
revealed from the photoemission, spin-lattice relaxation and neutron scat- 
tering experiments (see, e.g., Ref. [2]). However, this approximation has 
two serious shortcomings: (i) deviations from the Neel state are presumed 
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to be small and (ii) the approximation violates the rotation symmetry of 
the paramagnetic state. 

In this paper we try to overcome these difficulties by using the descrip- 
tion in terms of spin and hole operators and the continued fraction repre- 
sentations for Green's functions. Self-energies are calculated by decouplings 
corrected with the use of the constraint of zero site magnetization in the 
paramagnetic state. The Hamiltonian of the t-J model reads 

H= *nmaL a m CT + ^ J! Jnm ( S n S m + Sil^m 1 ) + ^^X n , (1) 

nmcr nm n 

where a ncr = |n<r)(nO| is the hole destruction operator, n and m label sites 
of the square lattice, a = ±1 is the spin, jno") and |n0) are site states 
corresponding to the absence and presence of a hole on the site. In these 
notations the spin-^ operators are written as s n = \ J2a c|ncr)(ncr| and 
s n = |ncr)(n, — a\. The hopping and exchange constants t nm and J nm are 
expected to be nonzero only for nearest neighbor sites for which they take 
the values t and J. fi is the chemical potential and X n = |nO)(nO|. 
The hole and spin Green's functions are determined as follows: 

G(kt) = -ie(t)({a ka (t),aij, D(kt) = -i0(t)([4(t), (2) 

where angular brackets denote canonical ensemble averaging, ako- = iV -1 / 2 
^2 n exp(— ikn)a no - and s k = iV -1 / 2 J2 n exp(— ikn)s n with the number of 
sites N. To calculate Green's functions (2) we use their continued fraction 
representation [3, 4], 

«A,l4»« = — ' ( 3 ) 

CO - En 



<jj — Ei 



where {{Aq\A)) 

w is the Fourier transformation of one of the above Green's 
functions and the operator An is equal to or s k . The coefficients Ei and 
F-i of continued fraction (3) are determined by the recursive procedure 

[A n , H] = E n A n + A n+ i + F n -iA n -i, 

(4) 

E n = |K, H] ■ Al\ \A n ■ A n \~\ F n _x = \A n ■ 4| \A n ^ ■ A^ 1 . 

The definition of the inner product \A ■ B' \ depends on the considered func- 
tion. In the case of the anticommutator Green's function G{kw) in Eq. (2) 
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this inner product is denned as ({A, B'}}, while for the commutator function 
D(kuj) it is ([A, B^]). However, in this latter case for the particular oper- 
ator Aq = s k the numerator of the continued fraction in Eq. (3) vanishes. 
Therefore the direct application of Eqs. (3) and (4) for the calculation of 
D{\tuj) is impossible. Instead we consider Kubo's relaxation function 

[•CO 

I4l4))t = oit) jf df<[tf(tV*k]>- (5) 

Continued fraction representation (3) with recursive procedure (4) is also 
valid for this relaxation function, if the inner product is defined as (A, B^) = 
i f£°dt exp(— rft){[A(t), #']), f] — > +0. On calculating the relaxation func- 
tion, Green's function is found from the relation 

= Kl-* k » + (6) 

For reasons of space we shall not go into details of the calculations 
and give only final results (a more detailed discussion will be published 
elsewhere). The spin Green's function reads 

n( , , 2a;(4,^ k )n(ka;)+4JC 1 ( 7k -l) 

D{kuj) = ^-2^n(M-^ ' (7) 

with the square of the frequency of magnetic excitations 

col = 16J 2 a|Ci|(l- 7k )(A + l + 7k ) (8) 
and the imaginary part of the polarization operator 

7j- roc 

ImII(ku;) = — J^/k'k / dkj'[riF(oj') — tif(v' — to)] 

LO k , J-oc 

xA(k' -k,u/ -w)i(k'w') (9) 

(the real part of II (ko;) is determined from the dispersion relations). In 
the above equations, (s k , sl k ) _1 = 4Ja(A + 1 + 7k ), A = C<il\C\\ + 
(1 - a)/(8a\d\) - 3/4, 7k = [cos(fc x ) + cos{k y )}/2, n F {u) = [exp(u/T) + 
T is the temperature, /k' k = 2iiV~ 1 / 2 ( 7k / - 7k /_ k )(s k , si k )~ 1 / 2 , and 
the hole spectral function A(ku) = — 7rImG(ko;). The parameters C p = 
iV _1 ^ k 7 k Ck with p = 1,2 and C k = J2n exp[ik(n — m)](s„ 1 s^ 1 1 ) ap- 
pear in the above formulas after the decoupling of correlations containing 
four spin operators. These correlations arise in the second step of recursive 
procedure (4) and the decoupling serves as the terminator for continued 
fraction (3). Following Ref. [5] we somewhat improve this approximation 
by multiplying the result of the decoupling with the vertex correction a. 
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The parameters C\, C2 and a are calculated from the above definition of 
C p and the constraint of zero site magnetization in the paramagnetic state, 

(4) = \ (i - x) - (s'ht 1 ) = 0, (10) 

which can be written in the form iV _1 Ck = 1/2 for small hole concen- 
trations x = (X n ). In this equation and in C p the value of Ck is determined 
by the equation 

C k = 4J|d|(l - TkK 1 coth[ Wk /(2T)], (11) 

which follows from Eq. (7). 

The parameter A describes a gap near the point (ir, ir) of the Brillouin 
zone in excitation spectrum (8). As follows from our calculations, the gap 
opens at x ~ 0.01 for T = and, in accord with the Mermin- Wagner theo- 
rem [6], for any nonzero T. The gap opening indicates the destruction of the 
long-range antiferromagnetic order and the establishment of the short-range 
order with the correlation length which is determined by the gap magni- 
tude. Notice that Eq. (9) is close in its form to the polarization operator 
obtained in Ref. [2] for the same model in the spin-wave approximation. 
However, the effective interaction constant /^'k^k 1 differs somewhat from 
the constant of the spin-wave approximation. 

As can be shown by analogous calculations of the transversal spin 
Green's function D ± (kt) = -id^is^^s^ 1 }, 

D ± (kt) = 2£>(ki), (12) 

as it has to be in the paramagnetic state. Thus, the approximations made 
do not violate the rotation symmetry of the solution. Equation (12) was 
used in the derivation of Eq. (11). 
The hole Green's function reads 

G(kcj) = 4>[w -£ k -n- E(kcj)]- 1 , (13) 

where <p = (1 + x)/2, the unrenormalized hole energy for moderate x 

e k w (l - 4a 2 Cf + 4Ci) ^H^, (14) 

and the imaginary part of the hole self-energy 

/oo 
duj'[nB(u)') + np(uj — a/)] 

xA(k-k',w- J)B (kV), (15) 
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Figure 1. The spin correlations Z(l) = 4|(s? (i0 nS( 00 n)| calculated for T/J = 0.5, 0.75 and 
1 in this work (open circles) and by the Monte Carlo method in Ref. [7] (filled circles). 
In both calculations a 32 x 32 lattice without holes was used. 



Figure 2. The hole spectral function AQ/llS) for the case of one hole in a 4x4 lattice, 
k = (7r/2,7r/2) and parameters J = 0.2t, r/ = O.lt. Our calculations for T = 0.02i and 
zero-temperature exact-diagonalization data of Ref. [8] are given in panels (a) and (b), 
respectively. 



with h kk , =32(^)- 1 t 2 {(7 k _ k; + 7k ) 2 + ( 7 2_ k; - 7 2)[(l +7k ,)/(l-7 k 0] 1 / 2 }, 
n,B(w) = [exp(cj/T) — and the spin spectral function B(kw) = — 7r _1 
ImD(ka;). It may be noticed that the general structure of Eq. (15) coin- 
cides with the respective equation in the spin-wave approximation [2] where 
the interaction constant is close to /i kk '. From Eq. (11) for an infinite lat- 
tice and small x and T we find d = -0.206734 and a = 1.70494 (this 
value of C\ = (sn 1 ^ 1 ) with n and m being nearest neighbors is close to 
the value obtained in Monte Carlo simulations [7]). With these parameters 
e k ss — 0.65t 7k . Thus the width of the unrenormalized band is much smaller 
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Figure 3. The hole spectral function along the (0, 0) — (tv, 0) direction of the Brillouin 
zone. The calculations were carried out on a 20x20 lattice for the parameters t = 0.5 eV, 
J = 0.1 eV, x 0.1 and T = 0. The Fermi energy is taken as zero of energy. In contrast 
to Fig. 2 the electron picture is used in this figure. 

than the bandwidth in the absence of correlations when the dispersion is 
described by 4tjk. This difference stems from the antiferromagnetic align- 
ment of spins which leads to spin flipping accompanying the hole motion. 

To check the validity of the approximations made we carried out a 
number of calculations which allowed comparison with exact diagonalisa- 
tion and Monte Carlo data. In Fig. 1 results of our calculations of the 
spin correlations are compared with Monte Carlo data of Ref. [7] . In Fig. 2 
the hole spectral function A(kuj) obtained from our calculations and by 
exact diagonalization in Ref. [8] are compared. In both calculations the ar- 
tificial broadening r] = O.lt was introduced to transform 5- functions into 
Lorentzians. As is seen from these figures, our calculations are in good 
agreement with the numerical experiments. 

The mentioned similarity of Eqs. (7), (9), (13) and (15) to the respective 
equations of the spin-wave approximation [2] at moderate doping leads to 
close resemblance in shape of the spectral functions obtained in these two 
approaches. As an example of the obtained results, the hole spectral func- 
tion is shown in Fig. 3 for values of t and J which correspond to parameters 
of cuprate perovskites. As seen from the figure, for moderate doping the 
hole spectrum contains two dispersive features for wave vectors near the 
boundary of the magnetic Brillouin zone. One of these features - a narrow 
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intensive peak slightly below (in the electron picture) the Fermi level u> = 
- is connected with the so-called spin-polaron band. In a lightly doped crys- 
tal the width of this band is of the order of the exchange constant J which 
is much smaller than the hopping constant t for parameters of cuprate per- 
ovskites. The spin-polaron bandwidth is characterized by the parameter of 
magnetic excitations because on the antiferromagnetic background the hole 
movement is accompanied by spin flipping. The short-range antiferromag- 
netic ordering is retained in moderately doped crystals. In these conditions 
a part of the spin-polaron band is preserved near the boundary of the mag- 
netic Brillouin zone. The second dispersive feature is a broad maximum 
the dispersion of which is characterized by the second energy parameter of 
the model t. The shape of this dispersion reproduces with some distortion 
the shape of the two-dimensional nearest-neighbor band. For wave vectors 
near the V point two broad maxima, one below and one above the Fermi 
level, are seen in Fig. 3. Both of them belong to the dispersion with the 
characteristic energy t. 

The Fermi surface derived from the spectral functions in the under- 
doped case is shown in Fig. 4a. The solid segments along the boundary of 
the magnetic Brillouin zone are formed by points where the spin-polaron 
band touches the Fermi level (notice that this band does not cross it). The 
Fermi level is crossed by the broader dispersive feature along the dashed 
curves. The spectral functions calculated for wave vectors on these curves 
are shown in Fig. 4b. On moving from (tt/2, 7r/2) to (vr/5, ir) on the bound- 
ary of the Brillouin zone the spin-polaron peak recedes from the Fermi 
level. The situation looks like no Fermi level crossing occurs near (tt/5,tt) 
and a gap opens between the energy band and the Fermi level. Only on 
another scale (compare Fig. 2) one can see that in this region the Fermi 
level crossing does exist. It is not seen in Fig. 4b because the maximum 
which crosses the Fermi level is completely lost at the foot of the more in- 
tensive spin-polaron peak and this peak is positioned somewhat below the 
Fermi level. Analogous behavior of the photoemission leading edge is ob- 
served in cuprate perovskites and is known as the photoemission pseudogap 
[9, 10]. Notice that the shape of the photoemission spectrum is determined 
by the spectral function and the correlation between their behaviors is ex- 
pected. In Fig. 4c the positions of the spin-polaron peak in the spectral 
functions and the leading edge of the photoemission spectrum measured 
[9] in underdoped Bi2Sr2CaCu20s+(5 are compared for wave vectors on the 
Fermi surface. The magnitude of the pseudogap calculated for x = 0.121 
and its symmetry are close to those observed experimentally. Moreover, at 
x = 0.171 which is close to the optimal doping concentration the calculated 
pseudogap disappears. This is also in agreement with experiment. 

Two distinct emissions were observed recently in the normal state pho- 
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Figure 4- (a) The Fermi surface for the underdoped case, (b) The hole spectral functions 
for wave vectors on the dashed curve of the Fermi surface in panel (a) . Curves from top to 
bottom correspond to k = (0.2tt, tt), (0.2tt, 0.9tt), (0.2tt, 0.8tt), (0.3tt, 0.7tt), (0.4tt, 0.6tt) 
and (0.57T, 0.5tt), respectively. T = 116 K, x = 0.121, t = 0.5 eV, and J = 0.1 eV. (c) The 
energy gap between the spin-polaron peak and the Fermi level along the dashed curve of 
the Fermi surface for x = 0.121 (squares) and x = 0.171 (diamonds) at T — 116 K . The 
position of the leading edge of the photoemission spectrum measured [9] in underdoped 
Bi2Sr 2 CaCu2 08 +( s along the similar Fermi surface is indicated by filled circles. 

toemission of the optimally doped Bi2Sr2- x La x Cu06+«5 in the vicinity of 
the point (ir,0) [11]. This crystal is a single CuC>2 layer material with suf- 
ficiently decoupled layers. Therefore the two emissions cannot be ascribed 
to the bilayer splitting. As mentioned, the hole spectral function of the 
two-dimensional t-J model has also two components and it is reasonable 
to compare them with experimental data. The photoemission spectrum in 
Fig. 5(b) was obtained from the spectral function by convoluting it with a 
Gaussian to model experimental resolution and cutting off with the Fermi 
distribution. In accord with experimental conditions of Ref. [11], the Gaus- 
sian width 20 meV and the temperature 45 K were selected. In spite of 
some differences in relative intensities and widths of the spectral features, 
the calculated spectrum has the same structure as the experimental spec- 
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Figure 5. (a) The normal state photoemission spectrum of the optimally doped 
Bi 2 Sr 2 - I La :l; Cu0 6+ «s (T c = 29 K) at T = 45 K and k = (0.62tt, 0) [11]. (b) The photoe- 
mission spectrum calculated from the spectral function of the t-J model for t = 0.5 eV, 
J = 0.1 eV, x « 0.1, T = 45 K and k = (0.67T, 0). (c) The dispersions of the two features 
in the photoemission spectrum of IfeS^-zLa^CuOs+i [11]. (d) The dispersions of the 
maxima in the spectral function in Fig. 3. 

trum shown in panel (a). As in experiment, the maximum with a larger 
binding energy has a considerably larger width and stronger dispersion in 
comparison with the maximum with a lower binding energy. The experi- 
mental and calculated dispersions of the two spectral features are compared 
in panels (c) and (d). As is seen from the figure, they are similar and the 
binding energies are of the same order of magnitude. Recently two disper- 
sive features were also resolved in the normal-state photoemission of the 
underdoped and optimally doped Bi2Sr2CaCu20g+5 with two CuC>2 planes 
in the unit cell [12]. This splitting was connected with a coupling between 
the adjacent Cu02 planes. Since energy distribution curves and values of 
the band splitting in that work are similar to those observed in the single 
Cu02 plane material Bi2Sr2_ x La a: Cu06+ < 5, another possible mechanism of 
the band splitting in Bi2Sr2CaCu20s + 5 - the above-discussed peculiarity 
of the spectrum of the two-dimensional strongly correlated system - has 
to be also taken into account. It counts in favour of this latter mechanism 
that the Fermi surface obtained in Ref. [12] (see Fig. lg there) is very close 
to the Fermi surface in Fig. 4a. 

In summary, the t-J model of the CuC>2 plane was investigated with 
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the method which retains the rotation symmetry of the paramagnetic state 
and accounts properly for the kinematic interaction. The observed unusual 
properties of cuprates - the photoemission pseudogap and the band split- 
ting - were shown to be inherent in the two-dimensional strongly correlated 
electron system. 
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